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Abstract We investigate the joint description of the interest-rate term stuctures of 
Italy and an AAA-rated European country by mean of a -here proposed- correlated 
CIR-like bivariate model where one of the state variables is interpreted as a benchmark 
risk- free rate and the other as a credit spread. The model is constructed by requiring the 
strict positivity of interest rates and the asymptotic decoupling of the joint distribution 
of the two state variables on a long time horizon. The second condition is met by 
imposing the reversibility of the process with respect to a product measure, the first 
is then implemented by using the tools of potential theory. It turns out that these 
conditions select a class of non-afBne models, out of which we choose one that is 
quadratic in the two state variables both in the drift and diffusion matrix. We perform 
a numerical analysis of the model by investigating a cross section of the term structures 
comparing the results with those obtained with an uncoupled bivariate CIR model. 

Keywords Interest rates ■ Multidimensional CIR processes ■ Potential theory 

JEL Classification E43 

Mathematics Subject Classification (2000) 62P05 • 60J45 



1 Introduction 

The difficulty to model the evolution of the term structure of interest rates is witnessed 
by the existence of a large number of models present in the academic literature and in 
the financial practice, see e.g. [3]|22) for a review. Broadly speaking, these models can 
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be grouped in financially oriented arbitrage models, whose main objective is pricing 
interest rate sensitive contracts and measuring risk associated with the time evolution 
of the term structure, and economically oriented models that are embedded in more 
complex market equilibrium models. Among equilibrium models that of Cox, IngersoU 
and Ross (hereafter CIR) is certainly one of the most attractive. This model, introduced 
in OS], is characterized by two main properties: mean-reversion to an asymptotic state 
and absence of negative interest rates. Moreover, as Gaussian- like models {i.e. models 
founded on Ornstein-Uhlenbeck processes) generally develop numerically relevant tails 
in region of negative interest rates with growing time horizons, the CIR formulation is 
particularly popular in financial applications having as underlying portfolios composed 
of government bonds and long time horizons, such as the strategic asset allocation 
of life insurance segregated funds. However, well known limits of the CIR model are 
that the term structure can assume (see, e.g. [1^) only the following three shapes: 
monotonically increasing, monotonically decreasing and humped {i.e. increasing to a 
maximum and then decreasing), the need to allow the model parameters to vary with 
time in order to capture the observed evolution (see, e.g. [3]), and the difficulty to 
describe simultaneously all types of interest rate sensitive contracts, such as interest 
rate swaps, caps and swaptions (see, e.g. [15]). Moreover a single factor model is unable 
to describe simultaneously the evolution of the term structure of real and nominal 
interest rates. 

AH the above difficulties lead quite naturally to multi-factor extensions of the basic 
univariate CIR model. For example, already in the original model proposed by Cox, 
Ingersoll and Ross in [6], the instantaneous nominal interest rate is a linear combination 
of two independent state variables, the real interest rate and the expected instantaneous 
infiation rate, each evolving in time according to univariate diffusion processes, thus 
realizing the stochastic version of the well-known Fisher equation. Another example 
is the two-factor extension proposed by Longstaff and Schwartz [18j , where the two 
factors are used to express the short rate and its volatility. A different interpretation 
proposed for the two factor model is that the factors are linked to the short and long 
(w.r.t. the maturity of the contract) rates, as in the Brennan and Schwartz model 
[2]. A three-factor extension has also been considered and empirically investigated, 
among others, by Chen and Scott U on U.S. market data. The three factor setting 
is often motivated by the findings of Litterman and Scheinkman [17] according to 
whom the empirical description of the intertemporal variation of the term structure 
needs the use of three factors: the general level of interest rates, the slope of the yield 
curve and its curvature, that is associated with the volatility. For the euro market, a 
recent empirical investigation of the term structure evolution [20] has shown that two 
factors are sufficient for a description of the data with mean squared error Gaussianly 
distributed with about 10 basis points dispersion around the observed values. 

The aim of this paper to investigate the interest-rate spread between the Govern- 
ment debt of two selected European Union member states, Germany and Italy, in the 
hypothesis that the spread refiects the different market opinions of their respective 
credit quality. Both countries, together with France, are known to possess the most 
liquid and high-volume Government bond markets in Europe which provide observa- 
tions for a broad maturity range so that it seems reasonable to assume that the impact 
of liquidity premia in bond prices can be safely neglected. To model the joint term 
structures of interest rates we introduce a two-factor CIR-like model where one of the 
factors is interpreted as a benchmark risk- free rate and the other is a credit spread. In 
this sense the model follows the fractional recovery approach of Duffie and Singleton 
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[101 §7.2], although -as discussed later- our model is not afhne in the state variables. 
In fact, since it is natural to expect that the same macroeconomic factors affects both 
the level of interest rates and credit spreads, it is unclear to what extent a two indepen- 
dent factor model could describe the joint behaviour of the Italian and German rates. 
This issue is particularly relevant in the measurement of risk measures on portfolios 
composed by Italian and AAA-rated {e.g. German) government bonds. 

The two-factor model investigated in this paper is costructed according to the 
following requirements. 

First of all, from the financial point of view, a fundamental requirement of nominal 
interest rates modelling is to avoid negative interest rates. In the univariate CIR model 
this is guaranteed by the choice of the stochastic differential equation. However, since 
interest rates are expected to be strictly positive, it is also important to establish under 
which conditions on the parameters the rates do not vanish. In the case of the single 
factor CIR model, this question has been solved, in a different context, by Feller [12j . 
obtaining a necessary and sufficient condition. More generally, the hitting conditions for 
one- dimensional diffusion processes have been completely characterized , see e.g. [21j . 
On the other hand, the multi-dimensional case is much less understood. For financially 
oriented multi-factor models this question has been partially addressed in [11]. From a 
mathematical point of view, potential theory methods, when applicable, are the natural 
tools to analyze the hitting conditions for diffusion processes [13]. In fact, they have 
extensively been used in several contexts, albeit - in our knowledge - not for financial 
applications. In this paper we analyze by these methods the hitting conditions for 
multivariate correlated CIR-like processes and apply the result in the costruction of 
the correlated two-factor model. 

The second important feature we require on the model is that the correlation 
asymptotically vanishes. More precisely, we impose that in the limit of infinitely far 
time horizon the joint distribution of the benchmark risk-free rate and the credit spread 
decouples into the product of two Gamma distributions. 

Finally, among the models meeting the above requirements, we select the "minimal' 
class, by further requiring that the drift and diffusion matrix are quadratic in the state 
variables. We shall refer to this model as the "asymptotically decoupling correlated" 
model or ADC model. 

As previously stated, we investigate the ability of the ADC model to capture market 
behavior by applying the model to the joint description of Italian and German term 
structures of interest rates, at a fixed calendar date. In the application the numerically 
demanding calibration of the ADC model, for which there is no closed form expression 
for discount factors, has been performed with the well known technique of simulated 
annealing, in the (fast) adaptive version developed by Ingber [14]. Finally, we compare 
the results obtained with the ADC model with those obtained with a "simple" bivariate 
CIR model with uncoupled state variables. 

Our finding is that, for the particular cross section here examined, the ADC model 
does not increase significantly the accuracy in the description of the two term structures 
with respect to the one achieved by the uncoupled bivariate CIR model. However, the 
predicted risk-neutral joint distribution of the two models are different. This suggests 
that -once risk premia are inferred from the analysis of time series- the "natural" 
distributions could be different, implying different values of risk measures for the same 
portfolios. 
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2 The univariate CIR process 

In this section we introduce the tools of potential theory by discussing in a self- 
contained way the condition for the strict positivity of the univariate CIR process. 
Fix a filtered probability space (12, T , .7-*, P) equipped with a standard Wiener process 
w. The univariate CIR process X = \^Xt t £ R+} is defined as the solution to the Ito 
equation 

dXt = K[e - Xt]dt + a ^/Yt dwt 

Xo =xo ^ ' 

where k, 8, a are strictly positive parameters and > is the initial condition. In the 
celebrated paper of Feller [12] it is shown that the transition probability density of Xt 
is given by 

pt(a;o,x)=ce-("+'') (^^y \,^i{2V^) (2) 

where 

jr, u — CXq e , V — CX, V = — (3) 



a2(l- 

and la is the modified Bessel function of the first kind of order a. Notice that at fixed 
time t the dependence of pt on x is only through d, while c, u and v are constant 
parameters. The parameter v controls the behaviour of the probability density in ((2)1 
as a;J.O: for u > 1 pt{xo,x) vanishes, for v = 1 it converges to ce~", while for < < 1 
it behaves as x'^~^ and therefore it has a integrable singularity. In fact, the properties 
of the modified Bessel function for q > — 1 are such that Ia(y) is real and positive for 
any y > and that in the limit y [Q one has, see e.g. [TJ pp. 374], 

where r{z) = t^~^e~*dt is the Euler gamma function. The cumulative distribution 
function F{x,t) = f'xQ[Xt < x] is obtained by integrating ^ 

Ft{xo,x) = J^ dypt{xo,y)^ j^^dye-^''+y^(^^ h^i{2^) = {2cx-2v,2u) 

(5) 

where x^(x;n,\) is the cumulative distribution function of a non-central chi-square 
distribution with n degrees of freedom and non-centrality parameter A. 

The long time behaviour of the real random variable Xt is given by the probability 
density ■Ki,{x) — lim(_,oo pt(a;o, a;). It is simple to check that this limit is independent 
on the initial condition xq and nu is just the density of a Gamma distribution with 
parameters u — 2k9/(j'^ and uj — u/6, namely 

TT^{x) = Nx^-^e-'^'' ^ Nx^-'^e-"''/^ 

.fA" 1 (6) 



We give now a potential theoretical proof of the classical result, again due to Feller 
[12j . that the CIR process hits the origin iff u > 1. Referring to [13] for an exhaustive 
treatment, we recall the basic notions of potential theory of reversible Markov process. 
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The generator of the process X, solution to ([T]), is given by the following operator 
defined on smooth functions on K+ such that /'(O) = (this condition corresponds to 
the Neumann boundary so that the origin is a reflecting barrier) 

Lf{x) = ^a^xf'ix) + k{9- x)f{x) (7) 

A straightforward computation shows that L is symmetric in L2{^+,dni,), where 
Tiu is the Gamma distribution given in ([S}. Note that we use the same notation for the 
Gamma distribution and its density. 

The generator L can be written in the explicit self-adjoint form as 

Lf{^)^^^AM^)xnx)\' (8) 
SO that the corresponding Dirichlet form is 

rOC 2 rOC 

D{f):=- dMx)fix)Lfix)^^ dMx)xf'{xf (9) 
Jo ^ Jo 

We also define the qudaratic form D\ by 

/■oo 

D,if) = Dif)+ d^,{x)f{xf (10) 
Jo 

By standard theory, see e.g. [13] , the form defined by (|9]) is closable and the associated 
Hunt process is the solution to ((T}. We shall denote by Di also the closure of the form 
defined above and let ©i be its domain. 

We now recall that the capacity of an open set O C R+ is defined as 

cap(O) - inf Diif) , To ~ {f G : f{x) > 1 , x e O} (11) 

For an arbitrary set B C R+ the capacity of B is finally defined as 

cap(B) := inf cap(O) (12) 

O open : ODB 

A classical result, see e.g. [13] 4.3] of the potential theory for diffusion processes is that 
set with null capacity are never reached; such sets are called polar. 

Proposition 1 (Unidimensional Feller condition) For the Dirichlet form ((9|, the 
origin, i.e. the set {0}, is polar if and only if v > 1. 

Proof . It is convenient to introduce the quadratic form Dc, with c > 

POO 

Dc{f)=D{f)+c / d-K{x)f{xf (13) 
Jo 

and let cap^, be the associated capacity. Of course a set is polar with respect to cap^, if 
and only if is polar with respect to cap = capj^. We shall compute the capacity of [0, e) 
for a convenient choice of c. 

The minimizer for the variational problem defining capf,([0,e)) solves the equation 

r Lf{x)~cf{x) =0 x£{e,oo) 

{,/(x-) = l x(^[Q,e\ (14) 

I / G I/2([0,oo),d7r^) 
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The differential equation in (|14p reads 

^a^ X fix) + K{e - x) fix) - cfix) ^ (15) 



xfix) + u(l-^)f'ix)-^fix)^Q (16) 

that, modulo a change of scale, is a confluent hypergeometric differential equation [l]. 
Instead of using confluent hypergeometric functions, it is simpler to perform the change 
of variable fix) = gix) /nuix). A straightforward calculation gives 

xg"ix) +[2-v + v^]g'ix) + g - ^]gix) = (17) 

We now take advantage of the arbitrariness of c by choosing c = iva^) j (20) = k; in 
this way the solution of (|17|l satisfying the appropriate boundary conditions is simply 
given by 

^^^^^ G(x):=/°°dy^ (18) 
G(e) y 

Hence 

. '1 2 



(19) 



As cap^({0}) = lim^j^o cap^([0, e)) it is now easy to check that the capacity of the 
origin is null when v > 1 since all terms vanish as e J. 0; for v = 1 the capacity of 
the origin is still null since G(e) diverges logarithmically. Finally for < 1, by the 
asymptotic expansion of G(£), see e.g. [1], 

m-j^^^jh^n^-^)^'-{:eT\ ..0 (20) 

where we used p8|) and Fiu) — iu — l)rii/ — 1). It is now simple to check that 
cap,({0}) = A.(l-i.). 

n 



3 Independent CIR processes 

In this section we extend the results of the previous one to the case of independent CIR 
processes. Let (fi, JT, JT^, P) be a filtered probability space equipped with a standard 
n-dimensional Wiener process w = iw^, ■ ■ ■ ,w") and consider the uncoupled system 
of Ito equations 



l,...,n (21) 



dXl = K,[9, - Xl]dt + a, J XI dwl 



— Xq 



,4 > (i = l,...,n) 
and set Ui := 2Ki9i/af . Of course Xl > for any t £ R-j- and any i = 1, • • • , n. 
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The generator of the n-dimensional process X — {Xi, . . . , Xn) is given by the 
following operator defined for any smooth functions on R" such that dxif{x) = if 
Xi = 0, 

n 

Lf{x) = J2 [-^'^hrdx^xJix) + n,{e, - x,)dxj{x)] (22) 

i=l 

Since the processes in \21\ are independent it follows that L is symmetric in 
L2{^^,diiu), where v_ :— {vi,...,Vn) and dni^ is the product of n Gamma distri- 
butions with parameters Ui :— 2KjSj/(7^. Its density w.r.t. the Lebesgue measure on 
R" is TTi^{x) = riiLi^f'i (^i)) where tti/- is as in ((Gjl. 

Similarly to the one dimensional case, we address the question of which condition 
the parameters should fulfil so that the n-dimensional process X = {X^, ■ • ■ , X") does 
not hit the origin, i.e. when X]"=i does not hit zero. From the one-dimensional 
result it follows immediately that > 0,i = 1, . . . ,n , i.e. the processes X does not 
hit the coordinate axes, iff i^-^ > 1, i = 1, . . . , n. However a less stringent condition is 
sufficient to ensure that X does not hit the origin, namely iff X^^j^ vi > 1. This result 
is proven below firstly by a comparison argument and successively by using capacity 
theory. 

Proposition 2 (n-dimensional Feller condition) The n-dimensional process X := 
(X , ■ • • X") hits the origin with positive probability if and only if :— X^iLl ^' ^ 

Proof. We first show that if f — > 1 then X^ does not hit zero P-a.s. Let 

and introduce n independent processes as the solution to the 

equation 

dY^ = [k,9, - kYi]dt -f a, Jvidwl ,_ ^ 



j = l,...,n (23) 

Since X* > and > 0, i = 1, ■ • ■ ,n, by a standard result on one dimensional Ito 
equations, see e.g. [211 Thm. IX. 3. 7], for each i = 1, • • • , n we have X^ > P-a.s. It 
is therefore enough to prove that Y^ does not hit zero. Let 



1 1 
Zt:=j:^Y; ^■■=J:^-0 (24) 

From Ito's formula we get 



^=l ^» 



zt = z+ Er=i ^ f'ds [^,9, - kYs'] + 4 /* ^ 

' Jo ~[ Jo 

= z + J ds^!^- kZ.^ + Mt 
where Mt is a martingale with quadratic variation 



Yi 



(25) 



(M)t = V ^ / dsYl= dsZs (26) 

~1 cr, Jo Jo 



We thus see that Z solves, in the sense of the associated martingale problem, the 
stochastic equation ((T)) with xq = z, k = H, 9 — v/(2k), and a = 1. From the result on 
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the one dimensional CIR process discussed in Section[2]we then get tliat if 2R, v j (2k) — 
> 1 then the process Z is P-a.s. strictly positive. 
To show that if ^ < 1 then X hits the origin with positive probability we argue 
in a similar way. Let k :— min^— i ... „ Kj and define Y as the solution to (|23|) with H 
replaced by k. Then < a.s., i = 1, ■ ■ ■ ,n. Moreover, letting Zt :— Y27=l ^ 

by the same computation as above, we get that Z solves fl} with xq = z, n = k, 
9 — i^/{2k), and a = 1. The result follows. □ 



The Dirichlet form corresponding to the generator L in (|22p is given 

Dif) = -f dn,{x) fix) Lf{x) = \Y.T'i j dn,{x) x, [dc^Jix)] ' (27) 

As in Section O given c > we also define the Euclidean norm Dc 

Dcif) = Dif) + c / d^d^) f{x)^ (28) 



We shall also denote by Dc the closure of the form defined above and let V be its 
domain. In the next result we prove that if ^'i ^ 1 then the capacity of the origin 

vanishes. In the next Section we show how it implies an analogous statement when the 
CIR processes are not anymore independent but they are constructed with suitable 
correlations. 

Proposition 3 //X^ILi J'i ^ 1 then the origin {0} is polar for the Dirichlet form (|22l) . 

Proof. For notation simplicity we consider only the two dimensional case, n = 2 and 
choose 01 — 02 ~ o-i = a2 = 1. For e > set 

Ae ■- {{xi,X2) eRi : xi+ X2>e] (29) 

we shall construct a function fe ■ ml -> R+ with /e = 1 on \ Ae such that if 
vi + 1^2 ^ i then, for a suitable c > (hence for all c > 0) we have 

lim Dcife) = (30) 
e— >0 

by the variational definition of the capacity this implies capj,({0}) — 0. 

We choose fe{xi,X2) = fee ( ^^'^^.^ xi + -^^^TI^^'^) where he : R+ ^ M+ wiU be 
chosen later. To estimate Dc(fe) we perform the linear change of variables 

Vl V2 

■X\^ ' X2 



Vi +1^2 Vl+ V2 

(31) 

V2 , Vx 

; x\ H \ ^1 



so that 



xi = xi{r,s) = ^ uii ^ 7T V2 s 

VI + VI + v^ 

I VI+V2 VI+V2 

X2 = X2(r,s) = ^ V2r ^ ^ tt vi s 

v( + v^ VI + v^ 



(32) 
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We then have 

Dcife) = JJ^ dxidx2 Tr^^{xi)TT^2{^2)\^^l[dxife{xi,X2)]'^ + X2[dx2fs{xi,X2)]'^ + cfe{xi,X2f'j 



vl + q 



{ , "2 ,2 ^iir,srx2{r,sr'' + "J xi{r,sr-'x2{r,sr} 
+cl^lL±i^iViiV2 r dre-^^'+^^^heirf [^\s xi[r,sY'-\2{r,sY^-^ 

(33) 

By the change of variable s = ry we have 
dsxi{r,sY^~^X2{r,sy~^ 



1/1 +77'^ — 1 

— r 



(34) 

as weU as 

{ , N2 ^i(^' s)"^^2(r, s^-' + "J xi{r, sr-'x2{r, s^} 

-r''^+''^C2 

(35) 

for suitable constants Ci, C2 > depending only on ^1,^2- Therefore 

Dc{fe)= ^V^^i^2 rdTr'''+''--^e~^'''+''-^^{C2r[h',{r)f + cCih,{rf} 

(36) 

and we conclude the proof by choosing he as in the one dimensional case with parameter 
vi + 1^2 for an appropriate c > 0. □ 



4 A class of bidimensional correlated processes 

While multi-dimensional independent CIR processes have been widely employed to 
describe systems with mean reverting characteristics, correlated CIR-like processes are 
less popular. A general setting for multi-factor mean-reverting processes where interest 
rates and credit spreads are affine in the state variables has been investigated in the 
works of Duffie and Singleton and Dai and Singleton [7], where e.g. the number of 
state variables is three and 

dXt = k(6) - Xt)dt + Ey/S{Xt)dWt (37) 

where 6 £ Ri]_, and k, E and S{t) are 3x3 matrices, out of which the first two are 
constant while 5 is diagonal and affine in the state variables, so that it is possible to 
mix Gaussian and CIR-like processes. 
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Clearly, the choice of the correlation structure depends on the properties of the 
system to be modelled. We introduce here a class of bivariate processes where the 
correlation is such that the invariant measure of the joint process is equal to that 
of two independent CIR processes. As a consequence, in the asymptotic state the 
two processes decouple. We shall refer to this model as the asymptotically decoupling 
correlated (ADC) model. 

This approach is different to the introduction of a correlation on the underlying 
Wiener processes in pif) . Indeed we perturb both the martingale part and the drift 
terms in pif) in such a way the decoupling holds ast^oo. On the other hand for finite 
times the corrections can still be relevant. We thus analyze the class of bidimensional 
processes of the type 

"^(xa) -MXi,X2)dt + B{XuX2)d(^]^ (38) 

where wi and 102 are independent Wiener processes and we restrict the choice of A and 
B to second order polynomials in Xi and X2. Specifically, given the CIR parameters 
K,i,8i,ai > 0, and the correlation parameters > and 7 £ [— £2, £2] with 



A = 



AA ^ f m{l + PiX2}iei ^ Xi] + i^2a2Xi[02 - X2]\ 

A2) ^2(1 + P2Xi)[e2 - X2] + >^iaiX2[ei - Xi] ) 

B{xi,X2)B{xi,X2)^ = S{xi,X2) 



with S{xi,X2) = 



2 X (39) 

(JlXl-\- eiXlX2 'yxiX2 

JX1X2 CT|a::2 + £22:i2;2. 



and f3i = a,; = ^ i — 1,2 

Thus while ei and £2 measure the degree of correlation in the diagonal terms, 7 mea- 
sures the asymmetry between the two factors. If £1 = £2 = 7 = the process reduces 
to the independent bidimensional CIR process. 
Since 

det(S') = 2:12:2 (s:ie20-l + 2:2£icr| + ctictI) + xlxl{eie2 - 7^) (40) 

under the model conditions £1 > 0, £1 > and £i£2 > 7^ the diffusion matrix is 
definite positive. Moreover the eigenvalues ei,2 and eigenvectors ei,2 of 5 are 

1 



ei.2 2 



{afxi + 61X1X2) + {(T2X2 + £22:12:2) ± 

1 (41) 

\J[{^l^l + £1x1x2) - {(^1x2 + e2XiX2)f + 4(7xj/)2 



ei,2 



'"'0-12:1 + £12:12:2) - (0-22:2 + £22;ia;2) T 



2'^xy V 



\J [(0-12:1 + £12:12:2) - (o|x2 + £22:12:2)]^ + 4(72;i/)2 j , 1 



(42) 
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while, out of the possible representations of the matrix B, the simplest is obtained by 
canonical Cholesky decomposition 



B 



1 / alxi+£iXiX2 'yxiX2 



^Jcrfxi + £1X1X2 



^det{S) 



(43) 



The model as defined above shows two relevant features: (a) the two random vari- 
ables Xi{t) and X2{t) decouple as t — ^ 00; and (b) the processes {Xi, X2) does not hit 
the origin if V1+V2 > 1, where as for the uncorrelated process Ui = 2ni6i/ a-fii = 1, 2). 

We now show that the process {Xi,X2) is reversible with respect to the probabil- 
ity measure Hux,v-2-: which, as previously, is the product of two independent Gamma 
distributions of parameters vi and V2- In other words the generator L of [Xi,X2) can 
be written as 

Lf{x) = ^ ^ V ■ {^{x)S{x)Vf{x)) (44) 



where V denotes the gradient and V- the divergence. The proof (|44|) amounts to a 
straightforward computation. The explicit form of the left hand term in (I44|) is 

L = [Aidxi + A2dx2\ + i [Siidxixi + 28128x1x2 + 8228x2x2] (45) 
while the right-hand term is 
11 



2 ■k{x) 



V- 'K{x)S{x)Vf{x) 



1 1 

2^x) 

1 
2 

1 
2 



-^ja^i [Tv{x){Siidx^ + Sx2dx-2)f{x)]+ dx2[^{x){S2idx, + 8228x2) 

Siidxi ln7r(a;) -I- dx^Sn + 8218x2 ln7r(x) + 8x2821] dxif{x) + 
8228x2 ln7r(a;) -f 8x2822 + 8128x1 \m\{x) + 8x1812] 8x2 f{x) + 



^ [8118x1x1 + 28128x1x2 + 8228x2x2] fix) 

(46) 

In the expression above the first term in squared brackets is equal to Ai 



2 \ / ^l ^ ^ \ / 2 \ / ^^2 ~ 1 V2^ 

{(Jixi -I- eia;ia;2)( -r) + (f^l +£2:2) + (7a;ia;2)( -5-) +72^1 

Xi t>i X2 f2 



£1 7 

= Kl{6i - Xi) + Ki^X2(6l - Xi) + K2^{02 - X2) = Ai 

(47) 

and similarly the second term is A2, which completes the proof. 

By using tools from potential theory we now show that if i^i + 1^2 > 1 then the 
ADC process does not hit the origin. The Dirichlet form of the two factor correlated 
CIR process is 

75corr(/) = i / dlX,i,^2{x)^f{x)-8{x)^j(x) (48) 

We then have 
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Proposition 4 Let fi + U2 > 1. Then the origin {0} is polar for the Dirichlet form 
(|48ll. 

Proof . Let Sq be the diffusion matrix of two independent CIR processes with param- 
eters vi and V2, namely 

So(^i,^2)= f'''^'^' 2 ) (49) 

\ 02X2 J 

Recalling that the diffusion matrix S{xi,X2) for the two factors correlated CIR process 
has been introduced in (|39p and that ei£2 > 7^, a simple computation shows that 
for any x £ we have Sq{x) < S{x). This means that for each u G we have 
V ■ {S ~ So)v > 0. This bound translates directly to a comparison of the associated 
Dirichlet form, i.e. 

A„dip(/) < Dcorrif) (50) 

where -Dindip denotes the Dirichlet form of the two-factor independent CIR processes 
((9]). The statement now follows from the variational characterization of the capacity, 
see (|lip . and Proposition |3] □ 



5 Application to interest rate modelling 

We have used the ADC model to investigate the interest-rate spread between the Gov- 
ernment debt of two selected European Union member states, Germany and Italy, 
in the hypothesis that the spread reflects the different market opinions of their re- 
spective credit quality. Data for the German and Italian term structures are deduced 
from the average bid-ask prices of zero coupon bonds and strips of coupon bonds 
quoted on the market on Oct. 31, 2006 (time to). The corresponding interest rates 
have then been interpolated using a natural cubic spline at thirty equally spaced val- 
ues of time to maturity r = 1,2, ...30 years to build the two term structures of 
interest rates r ic(fo,to + t) {c ~ D,I) and the term structure of the spread 
r i-> s{to,to + t) = ii{to,to + t) ~ iu{to,to + r). 

The result of this procedure is reported in Fig. [T] where the German and Italian 
term structures are shown together with the term structure of the spread upwardly 
shifted by 3.65%. In this way it easier to compare the dependence on time to matu- 
rity of the three curves. In addition, the plot also shows the zero coupon swap term 
structure r >—> izcs{tQ,to + r) extracted with the standard bootstrap technique from 
the values of annual interest rates swaps (the swap rates used here are those versus 
the 6 months Euribor, computed using the 30/360 convention). Noticeably the spread 
between the zero coupon swap and the German curves is about 22 basis points and 
is fairly independent from time to maturity. On the contrary, the spread between the 
Italian and the German curves increases with time to maturity at a rate very similar to 
the German term structure. Bid-ask spreads on the term structures are not reported 
on Fig. [T] since they are all smaller than 3 basis points. 

The analysis is based on the following main assumptions: (a) there is no credit risk 
loading on German bond prices, noticeably rated Aaa by all main credit agencies; (b) 
the German term structure and the spread between the Italian and German rates can 
be described by two state variables, respectively the benchmark risk-free rate rt and 
credit spread st evolving in time as 

dX = d h) = A{rt, st)dt + Birt, st) d (""A (51) 
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where wi and wt are standard independent Brownian motions and the matrices A{rt, st) 
and B{rt,st) are defined within either 

1. a bivariate CfR model where the two factors are independent (model 1) so that 

\Ks(as — st)J \ agSt J 

2. the ADC model introduced in the previous section (model 2) for which 



A{rt,st) = 



Kr{l + l3rSt)[0r - Tf] + hisasrt[9s - St] 
Ks{l + Psrt)[Os — St] + KrarSt[9r — n] 



,2. 



D/ \D/ \T a/ \ f (^rrt + £r rt St J rt St \ (ros 
B{rt,st)B{rt,st) ^S{rt,st)=( 2 , ^°'^> 

with Pi — -4t Q-i — — i — r,s; 

a particular case of model 2 is that obtained for er = es = 7 = 0, when it collapses 
to model 1; we should refer to this particular case as the degenerate ADC model 
and use it for calibration purposes. 

Finally, since in this work the analysis has been restricted to a single calendar date, we 
assume that (c) the equations of the two models given above are expressed according 
to the risk-neutral probability measure, so that for the moment being we do not need 
to further specify the market price of risk. 

The formal setting is inspired to the well-known fractional recovery of market value 
setting of DufRe and Singleton [9], that in turn is inspired to the recovery rules of 
over the counter derivatives In this setting the prices at time to, P£i{X,tQ,T) and 
Pj{X,tQ,T) , of a German risk-free and an Italian risky zero coupon bond paying one 
euro in T are obtained by discounting at the risk- free rate r{t) and at the effective rate 
r(t) -I- s{t) (without loss of generality we have absorbed the fractional recovery rate in 
the definition of s{t)), that is 

Pd (X, t , T) = eQ [e- I :Ft] 

(54) 

Pl{X,t,T) = eQ [e-/*''['-(«)+^(")l*" I jT^] 

The prices can also be obtained by the hedging argument and according to the Feynman- 
Kac formula, by solving the partial differential equation 



dtPciX, t,T) + J2 MX, t)d^^PdX, t,T) + ^Yl S^J^^^^, Pc{X, t, T) = 

= [rit) + lcs{t)]Pc{X,t,T) 



,Pc(x,r,T) = 1 

with c = D,I and Ic 



(55) 



if c = _D 

1 if c = / 
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We recall that the one-factor CIR model admits an analytic solution for the price 
of the unitary zero coupon bond [6] 

P{Tt,t, T) = [e- '^("^'^"1 Tt] = fit, T) e-9(*'^) (56) 

where 



ed{T-t) 



depend on = 2kQ ja^ and on the so-called Brown-Dybvig parameters d and 0: 

d= \/ac2 +2cr2, i(d-f k) (58) 

By the independence of the two factors, model 1 admits an analytic solution for the 
price of the unitary zero coupon bonds 

Pl{X,t,T)^P{rt,t,T)P{st,t,T) ^ > 

For model 2 we have chosen to compute the expectation integrals in (I54p using the 
Euler-Maruyama scheme for the evolution of X and the Simpson quadrature rule for 
the (stochastic) discount factor. The time step h — 0.004 years and the number of sim- 
ulations A'' = 5000 have been chosen by requiring the difference between the numerical 
result and the analytic expression to be smaller than few basis points in the case of 
the degenerated correlated model. 

The two models have then been calibrated to the observed term structures. The 
calibration of model 1 is done in two steps: first the four parameters of the risk-free 
curve are determined on the German data, and then the four parameters describing the 
evolution of the spread are calibrated on the Italian curve, having fixed the risk-free 
ones. For the second model we have fitted simultaneously the two curves by minimising 
the sum of the squared differences between the values of the risk-free rates and the 
values of the spreads. The minimization is performed using the MatLab [19] fmincon 
routine for model 1, while for model 2 we have implemented (in C) a procedure using 
the fast adaptive simulated annealing algorithm of Ingber [T3] in order to speed up the 
computation by avoiding the use of time-expensive numerical derivatives. The results 
of the fits are reported in Table 1, while the differences between the fitted curves 
and the observed ones are reported in Fig. [2] Notice that the values of the correlation 
parameters Sr, £s and 7 are different from zero and the Vr.s parameters are both greater 
than one. 

The accuracy in the description of the German term structure is approximately five 
basis points, which has to be compared with the maximum bid-ask spread of about 
three basis points. On the other hand the description of the Italian term structure is less 
accurate with deviations ranging up to approximately twenty basis points. Although 
the two fitted parameter sets are different (for example 6r in model 1 is 5.44% while in 
model 2 is 4.55%), the two models show a very similar degree of accuracy. Analogously, 
the structure of the deviations shown in Fig. [2] is very similar, possibly indicating the 
presence of a missing extra factor to be included in the models. 

Figures 3 and 4 show the comparison between the distributions of rt and st at r = 5 
years and r — 30 years computed with the two models, both when the ADC model is 
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Table 1 Results of the fit for model 1 (bivariate CIR) and model 2 (ADC model); for cal- 
ibration purposes the parameters of the degenerated ADC model (second column) are fixed 
to those of model 1. In the second part of the table the values of i'r,s = 2^^ „ 6^ s/cr^ ^ and 
<^r,s = ^'r,s/Sr,s are reported. 





model 1 




model 2 


(bivariate CIR) 




degenerated 


non-degenerated 




3.46% 


^•o 


3.46% 


3.39% 




0.0398 




0.0398 


0.0636 


Or 


5.44% 


O-r 


5.44% 


4.55% 


(Tr 


4.55% 


Of 


4.55% 


3.87% 








(0) 


258 (0.3859) 


SO 


0.04% 


SO 


0.04% 


0.19% 


Ks 


4.0049 


Kg 


4.0049 


3.3345 


Os 


0.29% 


ds 


0.29% 


0.26% 


as 


2.58% 


as 


2.58% 


4.23% 






Mes) 


(0) 


114 (0.2046) 






7 





0.2800 


Ur 


2.0857 


Ur 


2.0857 


3.8728 


Us 


35.0593 


Us 


35.0593 


9.6116 


UJr 


2.608% 


UJr 


2.608% 


1.174% 


UJb 


0.008% 


Uls 


0.008% 


0.027% 



degenerate (showing the quality of the caUbration) and when the ADC model is non- 
degenerated. The plots show that in both cases the spread st has essentially reached 
the asymptotic distribution already at t = 5 years. On the contrary the risk-free rate 
rt shows a slower convergence, particularly in the case of model 2. The similarity of 
the deviations in Fig. [2] is presumably due to the "fast" decoupling of the two factors 
in model 2. 

As a final comment, we notice that although our results show that there is no 
substantial gain in the description of the term structures by using model 2 with respect 
to model 1, the inclusion of the correlations in model 2 modifies the asymptotic state. 
This is better appreciated in Fig. 5 where the joint density at r = 30 years of the 
two state variables is shown for model 2 and in Fig. 6 where the difference between 
the joint density of the two models is reported. Qualitatively, this is in agreement 
with the results found in analysing cross sections of a single term structure with an 
unidimensional CIR model, where it is well known (see, e.g., [3, p. 100) that different 
sets of parameters can provide very similar quality of description. To analyse to what 
extent the value of risk measures for portfolios composed of Italian and German bonds 
are affected, it is then necessary to specify risk premia for the ADC model and calibrate 
their values by the analysis of historical time series of bond prices. 
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Fig. 1 Term structure of the German (continuous line) and Italian (dashed line) interest 
rates. The two curves are obtained by interpolating the values derived from the bid-ask average 
quotations on Oct. 31, 2006 of zero coupon and strips of government coupon bonds using a 
natural cubic spline. The spread (dot-dashed line) between the two curves, shifted by 3.65% to 
ease the comparison with the German interest rate term structure, and the zero coupon swap 
curve (dotted line) are also plotted. 
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Fig. 2 Difference (in basis points) between the fitted term structure and the observed one for 
the bivariate model (continuous line) and the ADC model (dashed line) for the Italian term 
structure (upper plot) and the German term structure (lower plot). 
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Fig. 3 Risk neutral distribution of rt (upper plots) and st (lower plots) at r = 5 years (left 
plots) and r = 30 years (right plots) for the model 1 (continuous line) and the degenerate ADC 
model (histogram) estimated with 100000 Monte Carlo simulations. The asymptotic value of 
the two distributions are also drawn (resp. dashed and dotted lines). 




Fig. 4 Risk neutral distribution of rt (upper plots) and st (lower plots) at t = 5 years (left 
plots) and t = 30 years (right plots) for the model 1 and the ADC model (histogram) estimated 
with 100000 Monte Carlo simulations. The asymptotic value of the two distributions are also 
drawn (resp. dashed and dotted lines). Notice that for both models the distribution of st 
essentially coincides with the corresponding asymptotic one already at t = 5 years. 
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Fig. 5 Joint probability density distribution of {Tt,st) at t = 30 years for model 2 (obtained 
with 100000 Monte Carlo simulations). 




Fig. 6 Diflferonce between model 2 (obtained with 100000 Monte Carlo simulations) and model 
1 in the joint probability density distribution of (rt, st) at r = 30 years. 



